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ABSTRACT 

We study the mid-infrared properties of galaxies in 30 massive galaxy clusters at 0.02<z<0.40, using 
panoramic Spitzer/MlPS 24/im and near-infrared data, including 27 new observations from the LoCuSS and 
ACCESS surveys. This is the largest sample of clusters to date with such high-quality and uniform mid-infrared 
data covering not only the cluster cores, but extending into the infall regions. We use these data to revisit the 
so-called Butcher-Oemler effect, measuring the fraction of massive infrared luminous galaxies {K<K*+1.5, 
Lir>5x1Q 10 Lq) within r2oo, finding a steady increase in the fraction with redshift from ^3% at z=0.02 to 
~10% by z=0.30, and an rms cluster-to-cluster scatter about this trend of 0.03. The best-fit redshift evolution 
model of the form /sf oc (1 +z) n has n = 5.7+jg, which is stronger redshift evolution than that of L* 1R in both 
clusters and the field. We find that, statistically, this excess is associated with galaxies found at large cluster- 
centric radii, specifically r 5 oo<r<r2oo, implying that the mid-infrared Butcher Oemler effect can be explained 
by a combination of both the global decline in star-formation in the universe since z ~ 1 and enhanced star 
formation in the infall regions of clusters at intermediate redshifts. This picture is supported by a simple infall 
model based on the Millennium Simulation semi-analytic galaxy catalogs, whereby star-formation in infalling 
galaxies is instantaneously quenched upon their first passage through the cluster, in that the observed radial 
trends of f$F trace those inferred from the simulations. The observed f$F however lie systematically above the 
predictions, suggesting an overall excess of star-formation, either due to triggering by environmental processes, 
or a gradual quenching. We also find that f SF does not depend on simple indicators of the dynamical state of 
clusters, including the offset between the brightest cluster galaxy and the peak of the X-ray emission. This 
is consistent with the picture described above in that most new star-formation in clusters occurs in the infall 
regions, and is thus not sensitive to the details of cluster-cluster mergers in the core regions. 
Subject headings: galaxies: active — galaxies: clusters: general — galaxies: evolution — galaxies: stellar 
content 



1. INTRODUCTION 

The evolution of galaxies in clusters since z — 1 is ex- 
pected to reflect both changes in the raw ingredients, i.e. 
the properties of galaxies that have fallen into clusters from 
the field in this time, and the physical processes that have 
acted on those galaxies after infall. Early evidence for cluster 
galax y evolution was presented by But cher & Oemlerl (119781 
11984 hereafter B084), who found that the fraction of cluster 
members bluer than the cluster red sequence (ft,), by at least 
Ab_v=0.2 mag in the B-V rest frame, increases from zero in 
the local universe to ^0.2 by z~0.4 implying a rapid evo- 
lution of the cluster population over the last 5 billion years. 
Empirically the star-forming spiral galaxies found by B084 at 
z~0. 2-0.4 are mostly replaced by SO gala xies in local clus- 
ters dDressler et alJfl997t iTreu et al.ll2003l) . A simple inter- 
pretation is that clusters accreted blue gas-rich star-forming 
spirals at z>0.5- 1 and that these galaxies have been trans- 
formed somehow into the passive SOs found in local clusters, 
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by having their gas reservoirs depleted by one or more phys- 
ical processes within clusters, including for example ram- 
pressure stripping, starvation or harassmen t (for reviews see 
iBoseUi & Gavazzill2006t [Haines et al.|[2007h . 

However, more recent studies of the so called Butcher- 
Oemler (BO) effect and the evolution of the morphology- 
density relation have suggested that at least some of the ob- 
served evolution in fb is due to selection biases. Firstly, B084 
selected galaxies on optical luminosity (My) rather than stel- 
lar mass (or Mk), and were therefore susceptible to biases 
arising from low-mass spiral galaxies - the optical luminosi- 
ties of such galaxies are boosted by starburst activity, and 
thus they increasingly enter the samp les at higher redshifts 
dDe Propris et alj 120031: iHolden et alj 120071) . Similarly, the 
use of a fixed Ab-v and My for blue galaxy selection fails 
to take account of the youth (and thus brightness and blue- 
ness) of stellar populations in galaxies at higher redshifts rel- 
ative tojhejrtower redshift cou nterparts. To cou nter this ef- 
fect, lAndreon etail (l2006l) and lLoh et all (120081) have advo- 
cated the use of differential fe-corrections to associate blue 
galaxies with the same spectral classes of galaxies at all red- 
shifts, and to take into account the expected luminosity evo- 
lution of galaxies when defining the My (or better still Mk) 
limits. The above studies all conclude that once these selec- 
tion biases have been eliminated the redshift evolution of fb 
is significantly reduced, suggesting that there has been little 
evolution in the cluster galaxy population since z—l- Sim- 
ilarly, when using mas s-selected samples CM>4x 10 10 Mq) 
IHolden et ail (120071) and lvan der Wei et al.l (120071) find that the 
morphological composition of clusters and the morphology- 
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density relation has remained largely unchanged since z~0.8, 
as opposed to luminosity-selected (My < My +1) samples 
where significant redshift evolution is s een (Smith etail200l 
iPostman et ai]|2005l:lD~esai et al.ll2007l) . 

Another more subtle bias arises from the selection of the 
clusters themselves. B084's cluster sample was a heteroge- 
neous mixture of clusters identified from photographic plates 
or by their association with radio galaxies, and hence favored 
the inclusion of more extreme clusters at high-redshifts, in 
particular those with higher blue fracti ons as they would be 
easi er to detect and identify. Indeed iNewberrv et al.l (1 19881) 
and lAndreon & Ettoril (1 19991) showed that the high-redshift 
clusters of the B084 sample are much more X-ray lumi- 
nous and have higher velocity-dispersions and central sur- 
face densities than their low-redshift counterparts, suggest- 
ing that these selection biases could mimic evolutionary ef- 
fects, and illustrating the need for well-defined samples of 
clusters spanning a large redshift range. More recent stud- 
ies of X-ray selected clusters again find that both the scat- 
ter in the blue galaxy fraction and the trend with redshift is 
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2006) relative to tho se measured from optically-selected clus - 
ter samples (BQ84. IPePropris et alJl2007t iGoto et al.ll2003l) . 
These latter samples contain a large number of poor clusters, 
and although the influence of X-ray luminosity (Lx), veloc- 
ity dis persion (a) or clu ster ric hness on fh is still much de- 
bated dMargoniner et alj|2001t IWake et all 120051: f Goto et al l 
120031: iPopesso et al.l 120071: lAguerri et al.ll2007l) . there seems 
a general consensus that for massive (o->600kms _1 ), X- 
ray luminous clust ers there is little s catter or evolu t ion in 
f b out to z~0 8 (ISmailetalJ 119981: [Homeier et al.l 120051: 
Pogg ianti et al.ll2006t lAguerri et al.ll2007l) . 

Finally, the scatter in the blue galaxy fraction among the 
intermediate-redshift (0.15<z<0.3) clusters is large; values 
for individual clusters lying in the range fb ~ 0-0.2, i.e. val- 
ues typical of clusters at z=Q and z=0.5 respectively. This has 
often been broadly attributed to the dynamical status of the 
clusters, with actively merging clusters showing higher frac- 
tions of blue or star-forming galaxies than those apparently 
relaxed dMiller & Owedl2003lHMiller et al.ll2005h . However, 
equally this could be due to field galaxy contamination, as 
the blue galaxy fractions are estimated by statistically sub- 
tracting foreground and background objects by comparison to 
control fields, a process that becomes increasingly uncertain 
at higher redshifts as the level of field contamination rises. 
However, in an analysis based on only spectroscopically con- 
firme d members of 60 cluster s at z < 0. 1 1 covered by the 2dF- 
GRS. lDe Propris et alj d2007l) obtained a similarly large scat- 
ter in the blue galaxy fraction, including some clusters with 
fb > 0.4, and also found no correlations between fb and other 
cluster properties including richness, substructure, concentra- 
tion. 

The sensitivity of ISO and more recently Spitzer at mid- 
infrared (MIR) wavelengths has opened up a new window 
for studying star-formation in galaxy clusters. Th e most 
straightforward interpretation dKennicutt et al.l [20071) of the 
24/im emission is that it traces the dust obscured star- 
formation, while the observed UV or H a emission traces 
the unobscured one dCalzetti et alj 120071) . MIR observa- 
tions of clusters have reveale d a population of dusty star- 
forming cluster galaxies (e.g. iFadda et al.l 120001: iDuc et alJ 
20021: iBiviano et al.l 12001 iGeach et alj|2006l: iMarcillac et alj 
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ISaintonge et alj|2008h . Indeed, these MIR-detected star form- 
ing galaxies often have optical colors consist ent with the pas- 
sively evolving ear ly-type cluster galaxies (Wol f et all 120051 ; 
lHaines et al.l 12008b and would therefore b e mis sed by the 
traditional BO studies. Sainton ge et al.l d2008l) therefore 
combined the dusty star-forming cluster members, defined 
here as having mid-infrared star-formation rates higher than 
5M yr _1 , with the traditionally selected blue cluster mem- 
bers - all found within R < IMpc of the cluster centers - to 
show that the optical and IR populations contribute roughly 
equally to the observed BO effect out to z — 0.8. However 
several studies have shown that the dusty galaxies contribute 
significantly more to the integrated cluster star formation rate 
than the optically selected blue galaxies. For example star 
formation rates derived from the [Oil] emission line are typi- 
cally ~ 1 0- 30 x lower than ra t es estimated from I R luminosi- 
ties (e.g. iMetcalfe et ai]|200H IGeach et al.ll2009l) . However, 
these previous studies have analyzed either single clusters or 
at most a handful of heterogeneously selected clusters over a 
large redshift range, with 1-2 clusters per redshift slice, such 
that no statistical analysis of trends with cluster properties or 
redshift has yet been done in the MIR. 

In this article we revisit the BO effect, taking ad- 
vantage of recently obtained panoramic Spitzer MIR and 
ground-based NIR imaging of 22 clusters at 0.15<z<0.3 
from the Local Cluster Substructure Survey (LoCuSS; 
http://www.sr.bham.ac.uk/locuss; see also jj2^1| for more 
details) plus comparable data for Coma and Abell 1367 
at z=0.023, five clusters from the Shapley supercluster at 
z=0.048 covered by ACCESS (A Complete CEnsus of Star- 
formation and nuclear activity in the Shapley supercluster), 
and CI 0024+1 7 at z=0.394. In con trast to the only previous 
MIR BO study dSaintonge et al.ll2008l) . our large sample of 30 
clusters in total allows us to measure the scatter in the f$F- 
redshift relation as a function of redshift, and to explore the 
relationship between the scatter and simple indicators of the 
dynamical state of the clusters. We also use the very wide 
field of view of our data to explore the radial dependence of 
the MIR BO effect out to cluster-centric radii of /?^2Mpc. 
In summary, we present a statistical analysis of 30 clusters, 
aiming to probe the balance between (i) the evolution of field 
galaxies that fall into clusters, and (ii) the physical processes 
(some of which may be related to cluster-cluster mergers) at 
play within the clusters, in shaping the population of actively 
star-forming cluster galaxies. 

In $2] we summarize the data used in this paper, and the 
photometric analysis. The main results are then presented in 
5(3] followed by a summary and discussion in $4] Throughout 
we assume £1^=0. 3, £!a=0.7 and Ho=70kms" 1 Mpc" 1 . 



2. DATA 

We have assembled a dataset on 30 clusters at 0.02<z<0.4 
with NIR imaging extending out to the infall regions and 
panoramic S/?/fz<?r/MIPS 24/im photometry covering the same 
regions. Observational details for the clusters are listed in Ta- 
ble 1 . New observations of 27 clusters at z < 0.3 are described 
in § jj2.1| & [2.2| below. The data and photometry for the one 
re maining clust e r at z = 0.4 (CI 0024) were de scribed in detail 
bv lMoranetaLl d2007l) and lGeach et all d2006l) : in this paper 
we make use of the CI 0024 master catalog 8 and refer readers 
to the relevant papers for further details. 

8 available from httD://www.astro. caltech.edu/~smm/clusters 
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Table 1 

Observational details of the cluster sample 



Cluster 


z 


MR" 




Ngals 


fsF 


name 






(Mpc) 


«a:*+i 


.5) 


Abell 1367 


0.022 


1 


0.83 


47 


0.068!™« 


Coma 


0.023 


1 


1.50 


126 


°- 029+ oo,7 


Abell 3556 


0.048 


2 


0.70 


16 




Abell 3558 


0.048 


2 


1.29 


181 


0.017+« : »j{ 


Abell 3562 


0.048 


2 


0.91 


61 


o.o52;o:o« 


SC 1327-313 


0.048 


2 


0.91 


128 


031+ 0028 


SC 1329-317 


0.048 


2 


0.76 


67 


0.051 


RXJ1720. 1+2638 


0.160 


3 


1.53 


183 


°- 048+ om6 


Abell 586 


0.171 


3 


1.15 


225 


°- 033+ oo,7 


Abell 1914 


0.171 


3 


1.56 


210 


°- 028+ oo!4 


Abell 2218 


0.174 


4 


1.26 


219 


°- 032+ 0016 


Abell 2345 


0.176 


3 


1.05 


118 


°- 021+ om5 


Abell 665 


0.182 


4 


1.38 


233 




Abell 1689 


0.182 


3 


1.50 


259 


°- 043+ oo!4 


Z1883 /ZwCl 0839.9+2937 


0.194 


3 


1.11 


94 


0.095+og^ 


Z1693 /ZwCl 0823.2+0425 


0.223 


3 


1.00 


110 


0.03 1 + <}<$ 


Abell 2219 


0.225 


4 


1.49 


363 


°- 043+ oon 


Abell 1763 


0.228 


3 


1.22 


262 


— u.Ulo 


Abell 2390 


0.230 


4 


1.50 


303 


0.073+«-<g7 


RXJ2 129.6+0005 


0.234 


3 


1.23 


238 




Z2089 / ZwCl 0857.9+2107 


0.235 


4 


1.02 


57 


0.089!°;™ 


Abell 1835 


0.252 


3 


1.59 


392 


063+° 022 


Z348 / ZwCl 0104.4+0048 


0.254 


5 


1.00 


131 


n 1 9Q+0.048 
u - lz5, -0.044 


Z7160 / ZwCl 1454.8+2233 


0.258 


3 


1.13 


114 


049+° 036 

-0.023 


Abell 1758S 


0.273 


3 


1.38 


189 


052+ 029 
-O.02J 


Abell 1758N 


0.279 


3 


1.16 


313 


053+ 021 


Abell 689 


0.279 


3 


1.20 


180 


| ,-,+0.039 
U - 1J/ -0.033 


Abell 697 


0.283 


3 


1.51 


286 


076+° 025 
u - u/ -0.022 


Abell 611 


0.288 


3 


1.37 


252 


n 126+° 025 


CI 0024+ 17 


0.394 


6 


1.00 


149 


n 1AR+0.043 
°- 148 -0 03R 



"So urce of near-infrared data: ( 1 ) GOLDmine database I Gav azzi et alj 
[2003b ; (2) WFCAM AT-band data from ACCESS; (3) WFCAM JK-band 
data from LoCuSS; (4) NEWFIRM 7^-band data from LoCuSS; (5) 
UKIDSS; (6) Treu et al. (2003) 



2.1. LoCuSS 

LoCuSS is a multi-wavelength survey of a morphologically 
unbiased sample of 100 X-ray luminous galaxy clusters at 
0.15<z<0.3 drawn from the ROSAT All Sky Survey clus - 
ter catalogs (lEbeling et al.ll998ll2000l:lBohringer et al.l2004l) . 
The overall aim is to constrain the cluster-to-cluster scatter 
in the observable properties (e.g. X-ray temperature, inte- 
grated SZ-effect T -parameter, star-formation rate, far-infrared 
galaxy luminosity function) of massive clusters at low red- 
shift, and to interpret these observables in the context of hier- 
archical assembly, aided by gravitational lens ing probes of the 
distrib ution of dark matter in the clusters (e.g. lSmith & Taylor! 
2008). These analyses will, for example, deliver new con- 
straints both on the normalization, shape and scatter of mass- 
observable scaling relations required for precision cluster cos- 
mology, and on the evolutionary pathways of gas-rich field 
galaxies into passive cluster early-type galaxies. Early re- 
sults on thes e two compl e menta r y aspects of the survey can 
be found in IZhang et alJ (120081) . iMarrone et all (l2009t) and 
lHaines et al.l d2009h . 

The first batch of 30 clusters in our survey benefits from a 
particularly ri ch dataset, includin g: Subaru/Suprime-Cam op- 
tical imaging (lOkabe et alJl2009l) . Spitzer/MTPS 24/im maps, 
GALEX near- and far-ultraviolet (NUV/FUV) imaging, and 



near-infrared (NIR; J,K) imaging from UKIRT/WFCAM and 
KPNO-4m/NEWFIRM. All of these data embrace at least a 
half-degree field of view centered on each cluster, and thus 
probe the clusters out to ~ 1 — 2 virial radii. We have also 
been awarded 500ksec on Herschel as an Open Time Key 
Programme to observe this sample at 100 and 160/im with 
PACS. These 30 clusters were selected from the parent sam- 
ple simply on the basis of being observable by Subaru on the 
nights allocated to us. In principle these 30 should therefore 
not suffer any gross biases towards one type of cluster or an- 
other (e.g. cool core cluster, merging cluster, etc.). In this 
paper we analyze the 22 clusters from the full sample of 30 
for which both NIR and MIR data are in-hand. 

2.1.1. LoCuSS Mid-IR Observations 

Each cluster was o bserved across a 2 5 / x25 / field of view 
at 24/zm with MIP S dRieke et alJl2004h on board the Spitzer 
Space Telescope 9 dWerner et al.l 120041) . consisting of a 5x5 
grid of MIPS pointings in fixed cluster or raster mode (PID: 
40872; PI: G.P Smith). At each grid point we performed 
two cycles of the small-field photometry observations with a 
frame time of 3s, producing a total per pixel exposure time 
of 90s. The central 5'x5' tile of some clusters had already 
been imaged by Guaranteed Time Observations program 83 to 
a much deeper depth (~3000s/pixel); these data were com- 
bined with our 24/im data where available to give complete 
coverage of the entire 25'x25' field centered on each cluster 
in the sample. 

The 24/im data were reduced and combined with the Data 
Analysis Tool (DAT ) developed by the MIPS instrument team 
dGordon et al.ll2005T) . A few additional processing steps were 
also applied as described in Ega mi et all ([2006). The data 
were resampled and mosaicked with half of the original in- 
strument pixel scale (l."245) to improve the spatial reso- 
lution. The 24/im mos aics were analyzed with SExtractor 
dBertin & Arnoutslll996l) : following SWIRE (Spitzer Wide- 
area I nfra-Red Extragalactic Legacy Survey; Lonsd ale et alj 
120031) we estimated the flux of objects within an aperture of 
diameter 21", and applied an aperture correction factor of 
1.29. The flux detection limits and completeness of each mo- 
saic were determined by individually inserting 500 simulated 
sources for a range of fluxes and determining their detection 
rate and recovered fluxes, using identical extraction proce- 
dures. The sources used in the simulations were formed by ex- 
tracting isolated, high signal-to-noise and unresolved sources 
from the mosaic itself. From these simulations, we estimate 
that the 90% completeness limits of our 24/im mosaics are 
typically 400/iJy, and for individual clusters in the range 300- 
500/iJy, the variation due to changes in the background cirrus 
level. 

2.1.2. LoCuSS Near-IR Observations 

The same 22 clusters were observed either with WF- 
CAM dCasali et alj|2007l) on the 3.8-m United Kingdom In- 
frared Telescope (UKIRT) 10 in March-November 2008 or 
with NEWFIRM on the 4.0-m Mayall telescope at Kitt Peak 1 1 
on 16-18 May 2008. The WFCAM data were obtained using 

9 This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA (contract 1407). 

10 UKIRT is operated by the Joint Astronomy Centre on behalf of the 
Science and Technology Facilities Council of the United Kingdom. 

1 1 Kitt Peak National Observatory, National Optical Astronomy Observa- 
tory, which is operated by the Association of Universities for Research in 
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FIG. 1. — J-K/K color-magnitude diagrams for galaxies within 1.5r5oo 
of clusters A 1689, A 2219 and A 2390. These three clusters are shown 
because their spectroscopic redshift catalogs are the most complete among 
the 22 clusters from LoCuSS. Filled symbols indicate spectroscopically- 
confirmed cluster members, with green (red) colors indicating those (not) 
having L/r>5 x 10 10 Lq based on their 24/jm fluxes. Blue and magenta open 
symbols indicate foreground and background galaxies with redshifts respec- 
tively, while crosses indicate those galaxies without redshift information. The 
pair of sloping lines indicate the upper and lower boundaries of the color- 
magnitude selection used to identify probable cluster members. The dashed 
line indicates our faint magnitude limit of K* + 1 .5. 



the same observing strategy as used by the UKIDSS Deep Ex- 
tragalactic Survey dLawrence et al.ll2007l) . covering 52'x52' 
to depths of K~19, J ~21, with exposure times of 640s, pixel 
size of 0.2" and FWHMs-0.8-1.1". The NEWFIRM data 
consist of dithered and stacked J- and TT-band images with 
exposure times of 1800s, and FWHM~1.0-1.5", covering 
the instrument field of view 27'x27' with a 0.4" pixel-scale. 
Each individual exposure was astrometrically and photomet- 
rically calibrated using 2MASS stars in the field, and stacked 
using IRAF, producing mosaics also reaching depths of A"~19, 
/~21. 

Following [Haines et al.l d2009l) . probable cluster members 
were identified from the JK photometry (with J-K colors 
determined in 2" diameter apertures), based on the empiri- 
cal observation that galaxies of a particular redshift lie along 
a single narrow J-K/K color-magnitude relation, as shown 
in Fig. Q] for A 1689, A 2219 and A 2390. This relation 
evolves redward m onotonically with redshift to z~0.5 (see 
lHaines et alj|2009h . The NIR colors of galaxies are relatively 
insensitive to star-formation history and dust extinction, with 
the J-K color varying by only ^O.lmag across the entire 
Hubble sequence, and hence only a single sequence is seen, 

Astronomy (AURA) under cooperative agreement with the National Science 
Foundation. 



unlike in the optical where separate red and blue sequences 
are visible. This is demonstrated in Fig. [TJ by the 24/im- 
detected cluster members (green symbols) showing the same 
J-K colors as those cluster members not detected at 24/tm, al- 
beit with more scatter. Indeed the 24/im sources often lie well 
above the sequence, presumably reddened by dust, necessitat- 
ing the use of a selection box that is asymmetric around the 
"red sequence" of cluster members. 

The NIR color cuts are particularly efficient at removing 
foreground galaxies, with >80% of K<K*+\.5 galaxies spec- 
troscopically confirmed as having Zsp-z c iuster< -0.05(1 +z) 
(shown by blue open symbols in Fig.Q} excluded by the lower 
J-K color cut. However there remains significant contam- 
ination by background galaxies in the range z c iusier<z<0-5 
(open magenta symbols), driven partly by the red envelope 
of our color-magnitude selection function in an attempt to 
minimize the loss of heavily reddened cluster members (see 
upper sloping line in each panel of Fig. [TJ. We correct statis- 
tically for this contamination, using galaxy counts satisfying 
the same color-magnitude selection in two control fields. We 
analyz ed the UKIDSS-DXS Lockman Hole and XMM-LSS 
fields dLawrence et al.ll2007l) . based on JK photometry cov- 
ering 2.1 deg 2 in total, obtained via the same observing strat- 
egy, reduction and calibration pipelines, as our WFCAM data. 
These fields also have publicly availab le 24/xm MIPS ph o- 
tometry from SWIRE DR2 catalogs (lLonsdale et al.ll2003l) 12 
which are complete to 450/dy. 

With the exception of the clusters highlighted in Fig. Q] 
the spectroscopic redshift information on galaxies in clusters 
from the LoCuSS sample is currently very sparse. For the 
purposes of this paper we therefore define galaxies as being 
cluster members if they satisfy the color-cuts discussed above. 
The fraction of star-forming galaxies in each cluster, f$p , cal- 
culated in §|3]from these cluster galaxy catalogs are therefore 
net of the statistical field subtraction described above, and the 
error bars include a term to account for the uncertainties in 
this subtraction. 

2.2. Low-redshift subsample 

For Coma and Abell 1367 we used archival 24/im 
Spitzer/MIPS data covering 2x2 deg 2 in the case of Coma 
(PID: 83, PI G. Rieke) and 30'x30' for Abell 1367 (PID: 25, 
PI G. Fazio). Both datasets were obtained in scan mode, and 
are complete to ~330/xJy, corresponding to L/s~1.4xlO 8 L 
dBai et al.l [2006). Note that the 24/im mosaics do not ex- 
tend out to the aperture of 1.5r5oo adopted in this paper, 
we therefore estimate the 24/im fluxes of galaxies that lie 
outside the observed 24/im field of view based on their 
IRAS 60/xm fluxes and the empir ical relation /6o//24~10 
dSoifer. Neugebauer & Houckll 1987b . The IR AS Faint Source 
Catal ogue 60/xm completeness limit is 0.6 Jy (Beichman et al. 
1988), corresponding to L/« ~ 4xlO lo L0 for Coma and 
A 1367, or just below our L/« selection limit used in |3] K- 
band photometry and redshifts for both of these clusters we re 
taken from the GOLDMine 13 database dGavazzi et al.ll2003l) . 

For the five clusters forming the core of the Shapley super- 
cluster (A 3556, A 3558, A 3562, SC 1329-313 and SC 1327- 
312) at z=0.048, we use unpublished Spitzer 24/im MIPS 
imaging (PID: 50510, PI: C. Haines), UKIRTAVFCAM K- 
band photometry and redshifts from the ACCESS survey 

12 The SWIRE DR2 catalogs are availab le here: 
http://swire.ipac.caltech.edu/swire/astronomers/data_access.html 
15 available at http://goldmine.mib.infn.it/ 
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all of which cover a ~2.5 x 1 .5 deg 2 field of view. The 
24/im mosaics are complete to ^400/iJy, corresponding to 
L IR ~\ x 10 9 L Q at the supercluster redshift. The WFCAM K- 
band data has exposure times 300s, pixel-scale 0.2" and is 
90% complete to £~17.5 (K*+5.5), while our redshift cover- 
age is >95% complete to K*+l .5. At these low redshifts many 
of the 24/im sources are resolved, and so for those galaxies 
with NIR diameters greater than our standard 21" MIR aper- 
ture, we used a series of circular apertures of diameter up to 2' 
(45" for Shapley), and the optimal diameter chosen to encircle 
all of the galaxy's NIR flux. 

2.3. Bolometric Infrared Luminosities q/24/im Sources 

The bolometric luminosity, L IR , of each 24/im-selected 
cluster galaxy is estima ted from its 24/im flux using a 
range of infrared SEDs dCharv & Elba j|200l: IChanialll2003l: 
iDale et all 120011: lLagache et all I2003L l2004h . as a function 
of redshift following iLe Floc'h et al.1 ([2005). We classify as 
star-forming those galaxies with Ljr>5x 10 10 L©, which cor- 
responds to the 24/im completeness limits for our five high- 
est redshift clusters (0.27<z<0.29). Note that for all the re- 
maining clusters at z~0.2, this L/« limit lies in the range 550- 
1500/iJy, and hence is always well above the corresponding 
completeness limit of the cluster. 

Assuming that the 24/im emission is due entirely to star- 
formation, this limit corresponds to a star formation rate 
(SFR) of ^Moyr- 1 (lKennicuttlll998h . The choice of in- 
frared SEDs introduces a spread of Lm/S24^ m ratios, due to 
the range of dust temperatures assumed. This spread varies 
systematically with redshift, and is of the order 0.1 -0.2dex 
over 0<z<0.3. These systematic uncertainties are incorpo- 
rated into the error bars quoted in $3] More accurate estimates 
of Ljr will be possible in future with the availability of Her- 
schel 100/im and 160/im photometry, which will allow us to 
fully model the reprocessed thermal emission from the dust. 

Some of the 24/tm emission could be due to dust- 
enshrouded active galactic nuclei (AGN), as a signifi- 
cant fraction of infrared galaxies are known to be pow- 
ered by a combination of AGN activity and star-formation. 
At the typical infrared luminosities of our sample (0.5- 
2xlO n L0), the primary power source of the 24/tm emis- 
sion has be en shown to be star-formatio n in >90% of the 
cases (e.g. iGoulding & Alexanderl 120091) . There is also 
a strong almost-linear correlation between the 24/im and 
extinction-corrected Paa fluxes over four decades in luminos- 
ity among local infrared galaxies, the latter emission shown 
from HST/NICMOS imaging to be in the form of nuclear star- 
formatio n rings or more extended e mission from spiral arms 
or disks (lAlonso-Herrero et al.ll2.Q06h . 

For our low-redshift cluster galaxy sample, we can con- 
strain the dominant power source behind the observed 24/im 
emission, using both optical spectroscopy and the mor- 
phology of the 24/im emission (exte nded/unresolved). For 
this purpose we use SDSS DR7 ( Abazaiian et alj [20091) 
for Coma and Abell 1367 members, or ISmith etalT (120071) 
for galaxies in the Shapley supercluster, to place 24/im- 
se lected galaxies on the emission-line diagnostic diagrams 
of iBaldwin. Phillips & Terlevichl 0981). We only consider it 
likely that an AGN as the dominant contributor to the infrared 
flux if a galaxy is both optically- classified as an AG N and the 
24/im emission is unresolved, as lBrand et al.l d2009l) found in 
a study of 0.15<z<0.30 galaxies with /24>0.3mJy that the 
bulk of those galaxies classified optically as AGN also had in- 
frared colors indicative of PAH emission from star-formation. 



We find that all 13 infrared-luminous galaxies in 
Coma/A1367 have either extended 24/im emission or optical 
emission-line ratios indicative of star-formation. In Shapley 
supercluster, we find just three of 29 infrared-luminous galax- 
ies are classified as an AGN. 

Considering also X-ray bright AGN, iMartini et all (120061) 
studied X-ray luminous (Lx>10 41 erg s -1 ) AGN in eight clus- 
ters at 0.06<z<0.31, obtai ning an AGN fractio n of ~ 5% 
among Mr<-20 galaxies. lEastman et a l. (2007) found evi- 
dence that the X-ray AGN (Lx > 1 42 erg s" 1 ) fraction increases 
by a factor 20 from z=0.2 to z=0.6. We consider what fraction 
of these X-ray luminous AGN are also strong MIR em itters. 
For Abell 1689, which is in the IMartini et akl {2006) sam- 
ple, only one of our infrared-bright sources is found to be an 
X-ray AGN, while for Abell 1758, we find just two of our 
82 24/im detection s to coincide with an X-ray point source 
(lHaines et al.ll2009l) . 

3. RESULTS 
3.1. The Mid-Infrared Butcher-Oemler effect 

For each cluster, we compute the fraction of star- 
forming cluster members (fsF) within 1.5r5oo (^1.0r2oo; 
ISanderson & Ponmanll2003l) . after selecting by rest-frame K- 
band magnitude (Mk<M^+1.5), and bolometric infrared lu- 
minosity (L/r>5 x 10 I0 L Q ). We exclude the brightest clus- 
ter galaxy (BCG) fr om each cluster, due to their unique star- 
formation histories dLin & Mohrl 12004) and the direct link 
between BCG activity and the presence of cooling flows 
within clusters (Ed gelll991l). We model the redshift evolu- 
tion of K* as a Bruzual & Chariot! (120031) stellar population 
formed at z=A with an exponentially-decaying SFR of time- 
scale 1 Gyr, norma lized to match th e present-day value of 
M£=-24.60±0.03 dJones et al.ll2006l) for field galaxies. We 
note that this value is slightly brighter tha n those observed fo r 
galaxy clusters, e.g. Mt=-24.58±0.40 dBalogh et al.ll2001l) . 
M*=-24.34±0.01 dLin. Mohr & Stanfordl2004l) (for a=\ . 1). 
but we prefer to refer to a "global" value for M* K . The mea- 
surements of rgoo for clusters at 0.15<z<0.3 are obtained 
from analysis of Chandra X-ray observations, both from the 
archive and our own Cycle 10 observations (PID: 10800565; 
PI: G.P Smith); the details of this analysis can be found in 
ISanderson et al.l d2009t) . Values o f r^pp for the remainin g clus- 
ters are taken from the literature (ISanderson et al.ll2006l) . 

The value of fsF for each cluster is listed in Table \T\ 
the quoted uncertainties in clude the binom ial error calcu- 
lated using the formulae in Gehrels (1986), and uncertain- 
ties on the statistical subtraction of field galaxies. We 
show fsF versus redshift - the MIR Butcher-Oemler effect 
- in Fig. |2] finding a steady increase in f$F with redshift, 
from (/s F )=0.035±0.023 at z<0.05, to (/ 5/ r)=0.053±0.027 
at 0.15<z<0.25, and (/ SF }=0.085±0.038 at 0.25<z<0.29, 
where the quoted uncertainties are the rms scatter around the 
means. This trend is consistent with that found for the original 
BO study, whereby f b increases from 0.053±0.054 at z<0.08 
to 0.098±0.062 over 0.17<z<0.28, although we caution that 
given the rather diverse counting radii, cluster and galaxy se- 
lection criteria used, this may be coincidental. 

To quantify the redshift evolution of fsF we fit the following 
relation to the individual data points at z < 0.3 shown in Fig. [2 
/sf=/o(1 +z) n - The best-fit exponent is n=5.1 J [ 2 l \; the corre- 
sponding curve and la confidence interval are over-plotted 
on Fig. Extrapolating the best-fit relation (shown as lighter 
shaded region) reveals that it is consistent with the value of 
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FIG. 2. — The mid-infrared Butcher-Oemler effect. The estimated fraction of Mjc<MJ+1.5 cluster members within 1.5r5oo having L;r>5x 1O 1O L0. The 
dashed line indicates the best-fitting evolutionary fit of the form f$p =/o(l +z)" to the z<0.3 clusters, and the shaded region indicates the Icr confidence region 
to the fit, with the lighter colors showing the extrapolation of the fit beyond ; 

indicate less X-ray luminous clusters. The open squares indicate the values of f$ F taken from the Saintonge et al. ( 2008) z>0.3 sample. 



f SF that we measure for CI 0024 at z = 0.395. Note that our 
measurement of fsF = 0.15±0.04 for this cluster is consistent 
with that obtained by iGeach et alJ (I2006I) . The rapid rate of 
evolution does not depend strongly on the choice of limit- 
ing radius; repeating the calculation using galaxies selected 
within rsoo we obtain n=5.0±2 3, however the values of f$F are 
in this case systematically lower than their counterparts within 
1.5r 500 by ~20%. 

Our use of NIR color cuts could exclude some star-forming 
cluster members, whose colors in Fig.Q]show a much greater 
spread than their passive counterparts, resulting in a fraction 
lying outside our color limits. To test this possibility, we 
redo the analysis, this time without any color cuts. The over- 
all trend and scatter are unchanged, although as the correc- 
tion for field contamination increases, so do the uncertain- 
ties both in fsF for individual clusters, and the overall mea- 
sured evolutionary trend, for which we obtain «=6.5^ o- The 
slight increase in (/ SF }=0.099±0.036 at 0.25<z<0.29 ob- 
tained when we remove the color cuts, suggests that a small 
fraction (^10%) of 24/im cluster members may be lost by our 
color cuts, but that this is not significant. 

The measured level of redshift evolution is also consistent 
wit h the MIR BO- s tudy o f eight clusters over 0.02<z<0.83 
by ISaintonge et all (120081) . at least for their z<0.4 clusters. 
The absence of a statistical sample of clusters at z>0 - 3 clus - 
ters in either this study or that of ISaintonge et al] ((2008), 
makes it difficult to interpret whether the strong redshift evo- 
lution might extend to higher redshift. However their two 
clusters at z~0.8 (shown as open squares in Fig. O are in- 
consistent with an extrapolation of our best-fit evolutionary 
model at 2a. This suggests that the MIR BO effect might sat- 
urate beyond z~0.5, although we note that they use a smaller 
fixed 1 Mpc counting radius to estimate their f$F ■ 

A possible caveat to this analysis is that our Ljr threshold 
for selecting star-forming galaxies lies within the exponential 
region of the infra-red luminosity function (LF), particularly 



at z~0 dBai et al.ll2006L estimate l; r =3x10 10 L q for Coma) 
so that small changes in the threshold could produce large 
changes in /sf ■ We are sensitive here to only the most actively 
star-forming galaxies, missing a significant fraction of the 
population of normal star-forming galaxies, which are known 
to form a well-defined sequen ce in the specific-SFR/stellar 
mass plane (Noeske et al. 2007). As a result we might expect 
our values of fgp to be systematically below the fi, values ob- 
tained for the same clusters. Ho wever, when aver aging over 
the six clusters in common with ISmail et all (119981) . we find 
no significant difference between our /sf values and their //, 
values. We also note that only two of the low redshift clusters 
(Coma and A 3558) have X-ray luminosities comparable with 
the clusters from LoCuSS at 0.15 < z < 0.3. Strictly speaking 
this analysis therefore does not compare like with like at low 
and intermediate redshifts. However, the results in this and 
subsequent sections are robust to the exclusion of all clusters 
at z < 0. 1 except for Coma and A 3358, in which case we now 
obtain fsF=T-^-2'l within 1.5r5oo- Indeed the primary driver 
of our observed rapid evolution in /sf is that produced within 
the LoCuSS sample, rather than the comparison to the low- 
redshift subsample. This mis-match is therefore not a major 
concern, however it is unfortunate that there does not exist a 
suitable comparable set of panoramic mid-infrared observa- 
tions of rich clusters at z<0.1. 

3.2. Mid-IR Butcher Oemler Effect and the Global Decline in 
Star-formation 

We now investigate whether the redshift evolution of /sf 
is caused by physical processes operating within clusters, 
or simply reflects the global decline in star formation since 
Z~l, or a combination of both. The redshift evolution of f$F 
discussed above is stronger than that of the gl obal star for- 
mation rate. For example, Zheng et al] d2007l) showed that 
the specific UV+IR-determined SFR has declined by a factor 
^10 since z~l, independent of galaxy mass, while the spe- 
cific SFRs of local galaxies are ~ 0.3dex lower than those 
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FIG. 3. — As in Fig. [2] but taking into account the evolution of the mid- 
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ing galaxies as star-forming. Symbols are as in Fig. [2] 
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FIG. 4. — As in FigTJ] but considering only those galaxies within rsoo- 
Symbols are as in Fig.[2] 



of 0.2 <z<0.4 galaxies of the same mass. iLe Floc'h et alj 
(2005) find strong evolution in the global infrared LF since 
Z ~ 1, characterized by Z^ g oc(l + z) m with m=3.2^1,, simila r 
to them = 3.0 ±0.3 obtained by lPerez-Gonzalez et all (12005b . 
iNoeske et all (120071) also find that the median SFR at a fixed 
stellar mass of the entire sequence of star-forming galaxies 
shifts downwards by a factor of three from z=0.98 to z=0.36, 
corresponding to an evolution of SFRoc(l +z) 2 9 . 

Studies of the MIR LFs of nearby clusters find the shape 
and LJ R to be consistent with that of field galaxies, while the 
evolution with redshift to z— 0.8 of L* g oc(l + z) 3 2±07 is also 
indistinguishable from the field dBai et al.ll2008l) . Similarly, 
iFinn et ail d2008l) find that the Ha luminosities of individual 
cluster galaxies have declined by a factor of up to ^10 since 
z~0.75, comparable to that of field galaxies over a similar 
redshift interval. 

To quantify the excess evolution of /sf over the global evo- 
lution of star formation we therefore repeat the analysis in 
jj3.ll this time allowing the luminosity threshold above which 
galaxies are counted in the numerator of f$F to scale with red- 
shift as follows: L /s >5xl0 10 [(l + z)/(l+z max )]*, with £=3.2 
chosen to be representative of the results discussed above, 
and z max =0.288, corresponding to the highest cluster redshift 
to which we fit the redshift evolution model. This is analo- 
gous to the use of differential ^-corrections to identify "blue" 
galaxies w ith the same spectral classes of galaxi es at all red- 
shifts (e.g. lAndreon et alj|2005 lLoh et al.ll2008l) . rather than 
the fixed A B _ y =0.2 criterion of B084. We re-plot the MIR 
BO effect using this selection function in Fig. [3] and find 
(as expected) that the redshift evolution has largely disap- 
peared: /s^=0.072±0.044 for the seven clusters at z<0.1, and 
/ SF =0.074±0.033 for the 22 clusters at 0.15<z<0.3. Again, 
we fit a model to the data of the form /sf(z) oc (1 +z) n , ob- 
taining n=1.22^[J2, as shown by the dashed curve and error 
envelope in Fig. [3] 

This level of evolution is consistent with the expected in- 
crease in the fraction of g alaxies accreted by m assive clusters 
within the previous Gyr. iBerrier et alJ (120091) show that al- 
though the accretion rate of galaxies into massive clusters has 
remained constant since z~l, half of the cluster galaxy pop- 



ulation has been accreted since z~0.4, and so the fraction of 
recently accreted galaxies should appear to double from z~0 
to z~0.4. This suggests that the "residual" redshift evolu- 
tion seen in Fig. [3] albeit at modest statistical significance, 
is attributable to the fate of galaxies after they have fallen 
into clusters. Galaxy populations that have recently fallen 
into clusters are most easily identified in photometric obser- 
vational data at large cluster centric radii, thus overcoming 
the projection of the three-dimensional distribution of galax- 
ies onto the sky. To gain a rough idea of the location of the 
galaxies responsible for the redshift evolution seen in Fig. [3] 
we modify further the selection function, this time restricting 
the range of cluster-centric radii to <rsoo, again using the "dif- 
ferential ^-correction" approach to Lir selection. The results 
of this modified selection are shown in Fig. [4] with the best-fit 
redshift evolution model again shown as the dashed curve; the 
best fit model has n=0.32*} gl We also fit a model with just a 
single parameter - a redshift-invariant value of fsF - obtaining 
fsF =0.056±0.004 with a reduced chi-squared value of 0.72, 
confirming that the scatter of the data around this value can be 
explained simply by the observational uncertainties, without 
recourse to either intrinsic cluster-to-cluster scatter or redshift 
evolution. 

In summary these results are consistent with the global de- 
cline in IR activity in field galaxies since z=l. As field galax- 
ies fall into clusters some of them suffer an increase in IR ac- 
tivity at >r5oo, presumably due to star formation induced by 
gentle processes in their local environment at large radii. The 
amplitude of this increase, as measured via f$F , evolves with 
redshift at a rate which is not significantly different from zero, 
but nonetheless consistent with the fractional rate at which 
galaxies are accreted into rich clusters. Nevertheless, there is 
still significant cluster-to-cluster scatter, and so now we turn 
to whether this variation in /sf might be due to some global 
property of the cluster. 

3.3. Correlations with cluster properties 

A number of pr e vious studies (e.g. IWake et ail 120051: 
iPopesso et all 120071: lAguerri et ail 120071) have looked for 
correlations between the X-ray luminosity of clusters 
and the fraction of star-forming galaxies (or equivalently 
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X-ray emission. 



blue/emission-line), with the aim of ascertaining if ram- 
pressure stripping and/or starvation could be the dominant 
mechanisms for quenching star-formation in infalling galax- 
ies, and driving the observed SF-density relation. 

In Fig. [5] we plot the fraction of star-forming cluster galax- 
ies within r5oo against the cluster X-ray luminosity in the 0.1- 
2.4 keV band tak en directly from the ROSAT All Sky Survey 
cluster catalogs (lEbeling etalJl!998l l2000t iBohringer etal] 
2004), using the same evolving L 1R cut as described in 
53.21 We exclude here Abell 689 as its X-ray luminos- 
ity is dominated by a BL Lac. As previously mentioned, 
clusters in the LoCuSS sample are more X-ray luminous 
(Lx>2x 10 44 ergs _1 ) than all of the other clusters except Coma 
and A 3558. Overall there is no app arent trend of fsF 
with X-ray luminosity, in agreement with lWake et alJ (2005), 
which might seem to rule out ICM-related processes. How- 
ever, this may simply be a saturation effect, in that ram- 
pressure stripping is effective at stripping the gas in all in- 
falling galaxies even for the lowest X-ray lu minosity clusters 
(Lx~3xl0 43 ergs _1 ) in our sample. Indeed, Poggianti et al. 
(2006) find a strong anti-correlation between the fraction 
of emission-line galaxies and cluster velocity dispersion for 
er^SSOkms" 1 , but for richer systems there are no systematic 
trends. This may also exp lain the observed ne gative trend be- 
tween fb and Lx seen by Pope sso et al.l (120071) . as their cluster 
sample extends to much poorer systems than ours or that of 
IWakeetalJ(l2005l) . 

As discussed in 5Q3 it has been suggested that the level 
of star-formation activity in a galaxy cluster correlates with 
its dynamical status, with merging clusters showing in- 
creas ed activity with respect to their undisturbed counter- 
parts dMiller & Owenl2003l:lMiller et alJ2005HMetcalfe et all 
120051) . No single measure exists in the literature that unam- 
biguously identifies a cluster as being "disturbed", i.e. under- 
going a merger. Nevertheless, several measures do appear to 
correlate with cluster dynamical status, such as the presence 
or absence of a cool core in the X-ray temperature profile, the 
cuspiness of the density profile, or the offset of the brightest 
cluster galaxy (BC G) from the peak of X-ray emission (e.g. 
Smith et al. 2005b). In these cases merging clusters are typi- 



cally identified with non-cool core clusters with flatter density 
profiles, and BCGs with large offsets from the X-ray peak. 

In Fig. [6] we therefore plot the fraction of star-forming 
galaxies within rsoo versus the projected physical distance be- 
tween the BCG and the peak of X-ray emission, for the Lo- 
CuSS subsample, the latter measurement being taken from 
Sanderson et al. (2009). Interestingly there is no obvious 
correlation between the BCG offset and f SF , or alternatively 
when comparing the cuspiness of the cluster density profiles 
with fsF- We also obtain similar results when considering 
galaxies within 1.5 rsoo. Among the LoCuSS sample, neither 
Abel l 1758 or Abell 1914, w hich are known merging clus- 
ters dOkabe & Umetsull2008h . have high f SF 's, while Abell 
611 which appears to be a regular, relaxed cluster is among 
the few "active" clusters with f S F=0.l35. Indeed, the most 
prominent feature in Fig.|6]is the absence of clusters in the up- 
per right quadrant of the plot, i.e. with large values of fsF and 
large offset between BCG and X-ray peak. On the face of it, 
this is counter to the qualitative expectation that merging clus- 
ters contain more numerous star-forming galaxies than non- 
merging clusters, as suggested from the radio observations of 
mergin g clusters (not among our sample) by iMiller & Owed 
(2003). This may suggest that at least within rsoo many of 
the galaxies have already been stripped of their gas when they 
were accreted into the progenitor (presumably already mas- 
sive) clusters, and so are unable to undergo any starburst phase 
triggered by the cluster merger. 

3.4. Radial population gradients and the infall model 

A comparison of Figs. [3]&[4]indicates that for many of the 
clusters, the fraction of mid-infrared sources is lower within 
rsoo than 1.5rsoo, in some cases by a factor two. This is 
suggestive of the well known morphology-density and SF 
versus density grad ients seen in clusters, both at low- and 
high-redshifts (e.g. iDressler et al I [19971 iBalogh et al] 1 2000 : 
Elling son et all 1200 U lTreuetal.1 120031; ISmith et al.l 12005 ; 
lHaines et al.l2007l) . It has also been noted that the blue galaxy 
fraction in clusters depends strongly on the radius within 
which the measurement is made, with fy, systematically in- 
creasing with cluster-centric radius (e.g. Ellingson et al.l200Tt 
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IWake et alJl2005h . 

In Fig.|7]we show the composite radial gradients in the frac- 
tion of mid-infrared luminous galaxies (as defined in j j3. 2\ 
for clusters in three redshift bins: the seven z<0.05 clus- 
ters (red squares and dashed-lines); the low redshift half of 
the LoCuSS sample (0.15<z<0.23; green circles and solid 
lines); and the high redshift half of the LoCuSS sample 
(0.23<z<0.30; blue triangles and dot-dashed lines). The er- 
ror bars indicate the uncertainty in the mean value of f SF 
for each sub-sample in each bin and not the cluster-to-cluster 
scatter. Each redshift bin shows a clear increase in /sf with 
projected cluster centric -radius, although within the uncer- 
tainties we find no strong evidence for redshift evolution in 
the radial trends plotted in Fig. [7] However outside r^, there 
appears to be a modest excess of galaxies in the highest red- 
shift sub-sample. This could be due to an increase in the frac- 
tional accretion rate of infalling galaxies with redshift, and/or 
an increase in the level of triggered/enhanced star-formation 
in galaxies as they encounter t he cluster environm ent for the 
first time, such as observed by iMoran et alJ (l2005h for galax- 
ies at t he virial radius of CI 0024+ 16. Similarly. iGallazzi et alJ 
(2009) find a significant population of galaxies undergo- 
ing highly-obscured star-formation, preferentially located in 
intermediate-density environme nts typical of t hose found near 
the cluster virial radius, while Fadda et al. (2008) find en- 
hanced star-formation activity at r>rsoo along two filaments 
feeding Ab ell 1763. These t rends are also qualitatively con- 
sistent with Ellings on et alJ ([2001 ), who found a steepening 
in the population gradients in clusters at 0.30<z<0.55 rela- 
tive to those at 0.18<z<0.30, and can explain the difference 
in the levels of evolution seen in Figs. [3] and |4] 

To test our infall interpretation of the radial trends seen 
in Figs. [3] IH & we examined galaxies falling into 20 
massive clusters (M^^IO^M ^) from the Millennium sim- 
ulation ( Spring e! et alJ 12005). These simulations cover a 
(500/r'Mpc) 3 volume, producing DM halo and galaxy cat- 
alogs based on th e semi-analytic models (GALFORM) of 
iBower et alJ d2006l) for which positions, peculiar velocities, 
absolute magnitudes and halo masses are all provided at 63 
snapshots to z=Q, allowing the orbit of each galaxy with re- 
spect to the cluster center to be followed. We select member 
galaxies from these twenty clusters that have K<K*+1.5, and 



lie within 5 Mpc of the cluster center, and identify which ones 
are infalling into the cluster for the first time. The fraction of 
infalling galaxies increases approximately linearly with pro- 
jected radius from close to zero at the cluster center, to ~25% 
at ^500 and ^55% at 2r5oo, until we reach ^4r5oo, beyond 
which we would not expect to find any galaxies that have 
passed once through the cluster (Mamon et al. 2004). We find 
no redshift dependence f or fi n f a in ng (r / r^w) within the simula- 
tions, at least over the redshift interval 0<z<0.3. 

We then consider a simple toy model in which the 
star-formation in these infalling galaxies is instantaneously 
quenched once they pass through the cluster core. We should 
not expect that all infalling galaxies are star-forming, a sig- 
nificant fraction will already be passive either due to pre- 
processing within groups or through internal mechanisms 
such as AGN feedback. We therefore use the galaxies from 
the UKIDSS-DXS fields in §|2j]to estimate the fraction of 
galaxies in the field that would satisfy our selection crite- 
ria, obtaining / JF =0.25±0.03 for the 0.15<z<0.23 clusters 
and / SF =0.29±0.04 for the 0.23<z<0.30 clusters, and use 
these values to re-normalize the radial profile of the Millen- 
nium galaxies at 4r5oo- For the low-redshift bin, we mea- 
sure fsF from the 561 Mk<M% + 1.5 galaxies (taken from 
the NYU-VAGC; iBlanton et alJl2005l) from the SDSS DR7 
having redshifts z<0.1 also having 24/im photometry from 
SWIRE, obtaining f S F=0. 19±0.02. Note that we are not using 
the GALFORM-produced SFRs to identify star-forming galax- 
ies, the only model-based parameter we consider is the K- 
band luminosity. 

The re-normalized curves are over-plotted on Fig.|7]as dot- 
ted lines whose color corresponds to its redshift bin. These 
show the same general trend for /sf to increase monotonically 
with cluster-centric radius. This consistency supports our in- 
terpretation of the radial trends in the data as arising primar- 
ily from the infall of star-forming galaxies from the field and 
suggests that a simple scenario where infalling star-forming 
galaxies are quenched once they pass through the cluster for 
the first time (for example via ram-pressure stripping), is valid 
at least as a first order approximation to model the evolution 
of the cluster galaxy population. The main caveat to this pic- 
ture is that the observed trends in f$p appear to lie systemat- 
ically above the predictions from the simple infall model for 
all redshift bins over 0.25<(r/r5oo)<1.5. One possible cause 
is projection effects, which can produce dramatic increases in 
the fsF for individual clusters, due to the presence of line- 
of-sight filaments and groups associated with the large-scale 
structure in which the cluster is embedded, at physical dis- 
tances 5-20 Mpc from the cluster, and which are not included 
in our simple infall model. Alternatively, this excess may indi- 
cate star-formation triggered by environmental processes, or 
a rather more gradual reduction in star-formation than the in- 
stantaneous shut-down in star-formation modeled here. Note 
that in our model, we make no distinction between those 
galaxies which pass through the cluster core on almost ra- 
dial orbits, and those which travel on more circular orbits, 
never passing within 0.5 rsoo- In future articles we will use the 
spectroscopic redshift information that is becoming available 
from our MMT/Hectospec redshift survey to examine in detail 
the properties and evolution of the MIR-bright and UV-bright 
cluster galaxy populations within the context of the infall sce- 
nario. 

4. SUMMARY AND CONCLUSIONS 
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We have presented a study of the mid-infrared proper- 
ties of galaxies in a representative sample of 30 massive 
galaxy clusters over the redshift range 0<z<0.4, taking 
advantage of panoramic Spitzer/MlPS 24/im and ground- 
based near-infrared observations from the LoCuSS and 
ACCESS surveys. We revisited the Butcher-Oemler ef- 
fect, using our infrared data both to reduce uncertainties 
on photometric selection of likely cluster members, and 
to overcome the strong dust-obscuration (~10-30x) pre- 
viously identified in comparative optical/infrared studies. 
We mainly considered fsp, the fraction of massive cluster 
galaxies (K<K*+1.5) within 1.5r 50 o with L IR >5x 1O 1O L . 
We find that (/sf) increases steadily with redshift from 
0.035±0.023 at z<0.05, to 0.053±0.027 at 0.15<z<0.25, 
and 0.085±0.038 at 0.25<z<0.29, where the quoted uncer- 
tainties are the rms scatter around the means. This trend is 
consiste nt with the trends in found by previ ous optical stud- 
ies (e.g. lSmail et alJll998l:rA.guerri et al.ll2007l) . however, the 
clusters-to-cluster scatter in f$f at fixed redshift is roughly 
half that seen in fy. The lower scatter in our results is likely 
due to a combination of our use of NIR data to select likely 
cluster galaxies and the different physics probed by optical 
and IR data. 

We fit a redshift evolution model of the form /sf oc( 1 +z) n to 
the observational data, obtaining a best-fit value of n=5.7±j \. 
This level of evolution exceeds that o f L* IR in both clus- 
ters and the field (e.g. iBai et al.l [2008). We therefore re- 
peated our analysis taking into account the cosmic decline 
in star-formation by modifying our selection function thus: 
L, R >5 x 10 10 [(1 +z)/(l +Zmax)]*, with £=3.2 c hosen to be rep- 
resent ative of the results suc h as those of iLe Floc'h et alJ 
d2005l) and lZheng etal] d2007l) , and z max =0.288, correspond- 
ing to the highest cluster redshift to which we fit the redshift 
evolution model. The best-fit redshift evolution for this modi- 
fied selection is: n=1.22±J jg. Indeed, if we restrict the galaxy 
samples to the central region of each cluster (r<r5oo) then 
the fraction of cluster galaxies identified as star-forming from 
their mid-infrared emission remains constant at ^5% over 
0<z<0.4. This suggests that redshift evolution of /sf can 
be interpreted as a consequence primarily of the rapid evo- 
lution in the SFRs of field galaxies over this period, which 
are accreted onto the clusters at a constant rate, before be- 
ing quenched by cluster-related processes. However the small 
residual redshift evolution seen at r>rsoo after removing the 
global decline in star formation suggests that some new star- 
formation is triggered in clusters at large radii, presumably 
due to gentle processes at play in the local group environ- 
ments within which galaxies arrive in the clusters. 

Globally, there is therefore little evolution in the cluster 
population itself, as the infall r ate of field galaxies is expected 
to remain constant out to z^l (Berrier et al. 2009), while the 
efficiency of the processes which quench star-formation in 
the recently accreted galaxies (e.g. ram-pressure stripping) 
should not evolve rapidly. This lack of evolution in the cluster 
galaxy population to z^0.5 and beyond, is consis tent with re- 
cent studies looking at the Butcher-Oemler effect dSmail et al. 
1998tlDe Propris et alj2003l:lHomeier et alj2005l:IWake et al. 
20051; lAndreon et al. 200^and the morphology-density rela- 
tion dHolden et al.l 12007). All of these studies indicate that 
much of the apparent evolution in earlier optical BO studies 
was due to the use of optical luminosities to select galaxies, 
such that the trends were due to low-mass galaxies under- 
going starbursts in the higher-redshift samples (see Fig. 2 of 
iHolden et al.ll2007h . or biases resulting from the cluster sam- 



ple itself. 

The view that the bulk of star-formation in clusters simply 
represents recently accreted field galaxies, is consistent with 
the observed constancy of the shape of the UV and IR LFs 
(both L* and a) from clusters to the fie ld in th e local Universe 
dCortese et all 120051: IBai et al.l l2006h . while iMercurio et ail 
(2006) find no variation in the optical LF of blue (i.e. star- 
forming) galaxies with environment. Moreover, the cluster IR 
and Ha LFs have been shown to evolve in the sam e manner 
as field LFs, declining by a factor ~10 since z~l (IBai et alJ 
2008: lFinnet"aT]|2008l . 

We have investigated the effects of cluster properties on the 
level of star-formation in their member galaxies, finding no 
apparent dependence of fsp on either X-ray luminosity or the 
dynamical status of the cluster. The absence of correlation 
with Lx may seem to rule out ICM-related processes as the 
main route by which star-formation is quenched in dense en- 
vironments, but this may simply be a saturation effect, in that 
ram-pressure stripping is effective in stripping the gas in all 
infalling galaxies even for the lowest X-ray luminosity clus- 
ters in our sample. Studies which examine much poorer sys- 
tems find strong anti-correlations between fsF and Lx or a v , 
which then flatten off for the er>550kms _I systems compa- 
rable to those which make up our sample. The lack of corre- 
lation between f$p and the dynamical status of the clusters 
is surprising and see ms to contradict previous studies (e.g. 
Mille /&" Owenll2003l) . This may suggest that at least within 
rsoo, galaxies have already been stripped of their gas, and so 
are unable to undergo any starburst phase triggered by the 
cluster merger. However at larger radii there may still be 
enhanced activity during certain phases of cluster mergers, 
as the infalling galaxies should still be gas-rich. Expanding 
on this, composite population gradients show a smooth in- 
crease in the fraction of star-forming galaxies from close to 
zero in the cluster cores to 7-13% by 2 7-500- Through com- 
parison with numerical simulations, we find these gradual 
trends are consistent with a simple model in which the star- 
forming galaxies are infalling into the cluster for the first time 
(usually on highly radial orbits), and then quenched somehow 
once they pass through the cluster core, for example via ram- 
pressure stripping. Within rsoo there is no apparent evolution 
in the radial population gradients, but beyond rsoo we find a 
possible excess of 24^m-bright galaxies in the highest red- 
shift bin (0.23<z<0.30), suggestive of either enhanced star- 
for mation in th e clust e r infall regions simila r to that found 
bvlMoranetal.1 d2005l) . iGallazzi et ail d2009l) or lFadda etal] 
(2008), or an increase in the fraction of infalli ng galaxies, 
comp arable to that expected from simulations dBerrier et al.1 
120091) . 

In the future we will further develop these results us- 
ing spectroscopic redshifts for the cluster galaxy populations 
from our ongoing M MT/Hecto spec survey, p lus weak-lensing 
dOkabe et alj|2009l) and X-ray dZhang et al.ll2008l) data to in- 
vestigate in detail the relationship between the star-forming 
galaxy populations in clusters and the dynamical state of the 
host clusters. 
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